Schwarz-preconditioned HMC algorithm for two-flavour lattice QCD
نویسنده
چکیده
The combination of a non-overlapping Schwarz preconditioner and the Hybrid Monte Carlo (HMC) algorithm is shown to yield an efficient simulation algorithm for two-flavour lattice QCD with Wilson quarks. Extensive tests are performed, on lattices of size up to 32×24, with lattice spacings a ≃ 0.08 fm and at bare current-quark masses as low as 21MeV.
منابع مشابه
Schwarz-preconditioned HMC algorithm for two-flavor lattice QCD
The combination of a non-overlapping Schwarz preconditioner and the Hybrid Monte Carlo (HMC) algorithm is shown to yield an efficient simulation algorithm for two-flavour lattice QCD with Wilson quarks. Extensive tests are performed, on lattices of size up to 32×24, with lattice spacings a ≃ 0.08 fm and at bare current-quark masses as low as 21MeV.
متن کاملHMC algorithm for two-flavour lattice QCD: Schwarz-preconditioning with a one-dimensional domain decomposition
We study a variant of the Schwarz-preconditioned HMC algorithm. In contrast to the original proposal of Lüscher, we apply the domain decomposition in one lattice direction only. This is sufficient to reduce the condition number of the fermion matrix restricted to the domains compared with the full fermion matrix. For the same linear extension of the domain, less links reside on the boundaries o...
متن کاملSpeeding up the Hybrid-Monte-Carlo algorithm for dynamical fermions
In large scale simulations of lattice QCD with two flavours of mass-degenerate Wilson fermions, the mass of the fermions is still too large compared with the up and the down quark masses. Therefore a delicate extrapolation of the data is needed. Unfortunately, the numerical effort required for the Hybrid-Monte-Carlo (HMC) [1] algorithm increases as the quark mass decreases for, at least, three ...
متن کاملInstabilities in Molecular Dynamics Integrators used in Hybrid Monte Carlo Simulations
We discuss an instability in the leapfrog integration algorithm, widely used in current Hybrid Monte Carlo (HMC) simulations of lattice QCD. We demonstrate the instability in the simple harmonic oscillator (SHO) system where it is manifest. We demonstrate the instability in HMC simulations of lattice QCD with dynamical Wilson-Clover fermions and discuss implications for future simulations of la...
متن کاملSolution of the Dirac equation in lattice QCD using a domain decomposition method
Efficient algorithms for the solution of partial differential equations on parallel computers are often based on domain decomposition methods. Schwarz preconditioners combined with standard Krylov space solvers are widely used in this context, and such a combination is shown here to perform very well in the case of the Wilson–Dirac equation in lattice QCD. In particular, with respect to even-od...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2004